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1. Introduction 



The cause for the late-time accelerated expansion of the universe remains one of the most 
compelling problems in modern physics. Many schemes have been proposed to explain this 
phenomenon. Although dark energy scenario jl| is the most popular one among them, ones 
have considered some models based on infra-red modifications to general relativity (GR), 
such as scalar-tensor theories H, f(R) gravity p| and braneworld models M. In general, the 
resulting field equations are fourth order because the Ricci scalar is constructed from the 
second order derivatives of the metric, and this feature may lead to pathologies. Recently, 
an alternative model based on modified teleparallel gravity receives considerable attention. 
See §, |, 0, |, |, [HI III], H, HI [IJ, 0, 0, H, EDI in detail - Instead of describing 
gravitational interaction with curvature of the background spacetime by employing the 
torsionless Levi-Civita connection, one can explore the opposite way and resort to the 
Weitzenbock connection that has no curvature, in this case torsion will independently do 



the job that curvature does in GR |2C]. There are some terms in the modified Friedmann 
equation that can be identified as the effective dark energy to give rise to the accelerated 
expansion of the late-time universe. This paradigm boasts the significant advantage that 
the field equations are second order, as we shall see in the next section. 

In order to discriminate different models, we need to break the degeneracy of back- 
ground expansion history, and work in first order perturbation to find more information 



concerning different models [21, 22, |2^, 24]. The matter density perturbation which char- 
acterizes the inhomogeneities of the universe comes to the rescue. While different models 
may have the same background behavior, their linear growths of matter perturbation can 
be quite different. See, for example, [25]. In this paper, we derive the evolution equation 
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of the linear matter density perturbation and find that it takes the same form as the coun- 
terpart in GR, except that the effective Newton's constant is rescaled by a term related to 
the first derivative of f(T). Note that when the Lagrangian of gravitational part is torsion 
scalar T, it is equivalent to the Einstein-Hilbert Lagrangian of GR up to a divergence, 
and hence all behaviors of this theory reduce to those in GR, including the local Lorentz 
symmetry. While for more general f(T) gravity, local Lorentz transformation fails as a 
symmetry of this theory p6| , and it is expected that this equation may also provide us 
some hint about Lorentz violation. 

The outline of this paper is as follows. In Sec. 2, we briefly review the theoretical 
structure of teleparallel gravity and how it explains cosmic acceleration. In Sec. 3, we 
present the first order equations based on metric perturbations and vierbein perturbations 
respectively. Moreover, we derive the governing equation for matter density perturbation, 
then solve it numerically and compare our result to the counterpart of GR. Finally in 
Sec. 4, we discuss and summarize our results. 

2. A brief review of f(T) cosmology 

In Riemann-Cartan spacetime, the curvature tensor and the torsion tensor coexist. On a 
manifold, one can define a large number of connections, which differ from each other up 
to a tensor quantity. The assumptions of torsion-free and metric compatibility lead to the 
Levi-Civita connection, and this is the one on which Einstein general relativity is based. 
However, we are free to choose other connections, for instance the Weitzenbock connection 
which is defined by 

= e x A d v e A , (2.1) 

where e^x^) is a set of orthonormal vectors, which form a noncoordinate basis for the 
tangent space at each point on this manifold, and e A (x li ) is the dual vectors. The torsion 
tensor is given by 

T\ t u = - r*, = e A A (0 M e/ " (2-2) 

one can find that the curvature tensor and the covariant derivatives of e^x^) with respect 
to this connection vanish globally, therefore ga{x^) are absolutely parallel vector fields, 
and this theory is dubbed teleparallel gravity |20|] . In this formalism, the fundamental 
dynamical object is the vierbein field e fl A (x), and the metric tensor is obtained by a 
byproduct 

g^u(x) = r] AB e fM A (x)e u B (x), (2.3) 

and then the Levi-Civita connection can be defined in a natural way. The difference between 
these two connections is described by the contorsion tensor which takes the form 

K^ p = -^(T^p-T^p-Tp^). (2.4) 

Finally, we can define a torsion scalar as follows 

T = S p ^T"p U , (2.5) 
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where 

V = \(^% + 5»T a \ - b v p T^ a ). (2.6) 

T is the simplest teleparallel Lagrangian, which differs from Einstein-Hilbert Lagrangian 
only up to a boundary term p7|] . 

Similar to f(R) gravity, we can write down the Lagrangian for the gravity as a function 
of T. The full action reads 

1 = i / dAxe • t T + f ^ + / d4xe • 

where e = det(e M " 4 ) = yj—g, k 2 = 8irG(G is the Newton's coupling constant), and C m 
stands for the matter Lagrangian. Performing variation in this action with respect to the 
vierbein yields the equations of motion 

-le a A [T + f] + eP A T^pS™[l + f T ] + e-%(ee? A + fr] 

+ e p A SrhTdyT = ye^T p a , (2.8) 

where the subscript T' denotes the derivative with respect to the torsion scalar and T p a is 
the energy-monmentum tensor. Comparing to the equation in f(R) gravity which is fourth 
order differential equation, this equation has an advantage of being second order. So it is 
much easier for us to analyze it. 

From now on, we focus on the a spatially flat Friedmann-Robtson- Walker (FRW) 
universe only filled with dust-like matter. The metric is given by 

ds 2 = dt 2 - a 2 (t)dx 2 . (2.9) 

In this case the torsion scalar is related to the Hubble parameter H = dlna/dt by 

T = -6H 2 . (2.10) 

The background equations of motion become || 

H 2 = ^p- f --2H 2 f T , (2.11) 
6 b 

■ _ 1 6H 2 + f + 12H 2 f T 

H ~ -ll + f T -12H*fTr' (2 " 12) 

where p is matter energy density. Evidently, the last two terms at the right-hand side of 
Eq. fl2.ll ) can be explained as the effective dark energy whose energy density is given by 

Pde = T^(-f + 2Tf T ), (2.13) 
and the corresponding equation of state is 

w= x | (f-T-2Tf T )(f T + 2Tf TT ) 
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From the above equation, the effective dark energy becomes an effective cosmological con- 
stant if 

f-T-2Tf T = 0, (2.15) 

or 

f T + 2Tf TT = 0, (2.16) 
for all T, The solution of Eq.([U5|) is 

f(T) = —T + civ 7 — T- (2.17) 

But now (1 + fx + 2T/tt) equals to zero as well. So this solution does not provide an 
effective cosmological constant. Switch to Eq.( |2.16| ), the solution is 

f(T) = c 2 V^T-K 2 A, (2.18) 



where A is a constant. Substituting the above solution into Eq.( 2.13j ), we find that the 



effective dark energy density is simplified to be A which is nothing but a cosmological 
constant and the term C2V—T does not contribute to the effective energy density at all. 
To summarize, f(T) can be taken as a cosmological constant only when f(T) is a constant. 
However it is quite trivial. 

3. Growth factor 

The early universe was made very nearly uniform by an inflationary state. The origin 
of structure in the universe is seeded by the small quantum fluctuations generated at the 
inflationary epoch. These small perturbations over time grew to become all of the structure 
we observe. Once the universe becomes matter dominated primeval density inhomogeneities 
(Sp/p ~ 10~ 5 ) are amplified by gravity and grow into the structure we see today.In this 
section, we investigate how the matter density perturbation grows in f{T) gravity. We 
keep terms up to the first order in the perturbed vierbein field. For the sake of simplicity, 
we will work in Newtonian gauge, which is valid for f(T) gravity theory because it still 
preserves the principle of general covariance. 



This section is divided into three subsections. Firstly, we follow the approaches in [18] 
and define all the scalar degrees of freedom in the perturbed metric in Newtonian gauge. We 
find that this ansatz is too naive and problematic. In the second subsection, we start with 
a general perturbed vierbein field which includes more degrees of freedom and derive the 
correct evolution equation for the matter density perturbation. In the last subsection, we 
consider a concrete f(T) model and compare the growth of matter over-dense perturbation 
in f(T) gravity with that in GR. 

3.1 A naive ansatz for the perturbed vierbein 

Up to the linear order, scalar perturbations should decouple with vector and tensor per- 
turbations. The perturbed FRW metric can be written by 

ds 2 = (1 + 2(j))dt 2 - a 2 (t)(l - 2iP)5 lj dx i dx j . (3.1) 
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A naive ansatz for the perturbed vierbein can be written by 



e » A -{~ o\(i-m m )' (3 ' 2) 

Accordingly, the perturbed energy-momentum tensor takes the form 

5T v = ( ~ 5p -a~ 2 {p + p)div\ 
^ \ (p + p)div 5j5p y 

where p is the energy density, p is the pressure and v is the velocity potential. From the 
action in fl2.7|), the first order perturbations are governed by 



l + 4> 



E° : i-bp = -^(1 + fr) - + H<f>)(l + f T - 12H 2 f TT ), (3.4) 

E\: ^{p + p)d i v = -{d i i) + H^){l + f T ) + 12HHf TT d i ^, (3.5) 

K 2 

E°i: —{p + p)d l v = -(d i 4> + Hd t <i))(l + fT-12H 2 f TT ), (3.6) 



E l j(i=j): 

K 2 

—Sp = -(36#V + mH 2 H(j) + 12H 3 (p + 36H 3 ifj + 36HHifj + l2H 2 ^)f TT 

+ (3H 2 (/) + H<j) + 2Hcf) + 3Hij) + $)(1 + fr) + UAH 3 Hf TTT (iP + H<p)(3.7) 
E l 3 (i+j) : ^-0 = 0, (3.8) 



which are the same as those in [12, 18], and we also use E^a to denote the equation 
obtained from variation of the action with respect to e^ A . Note that d % = d^dj and 
d 2 = d' l di throughout this work. Comparing Eq.([3,5|) with Eq. (|3.6| ), we approach to an 
extra scale-independent constraint on (j), namely 

Hdii; = H(d i 7p + Hd i( j>), (3.9) 

if Stt 7^ 0. Or equivalently, there are the same number of degrees of freedom as that in 
GR, but f(T) theory leads to one more equation. It may lead to inconsistency. Besides, 
the scale-independent evolution of <j) in the above equation is incompatible with the the 
integrated Sachs- Wolfe effect. 

One may consider /tt = 0. However, if so, f(T) oc T and f(T) gravity is effectively 
reduced to GR. In order to solve this puzzle, we go to a more general ansatz in the next 
subsection. 

3.2 A general ansatz on perturbed vierbein in f(T) gravity 

Since the local Lorentz symmetry is broken down in f(T) gravity, extra degrees of freedom 
compared to GR should appear. In this subsection, we include all scalar degrees of freedom 
in the vierbein. Following [|29|], the perturbed vierbein is expressed in terms of unperturbed 
vierbein e^ A and a first-order quantity xa B as 

e/ = (^ + XBV- (3.10) 
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with e$ A = 5q and ei A = aSf. In this subsection, we only take into account the scalar 
degrees of freedom, which are encoded in xb A as follows 

XAB = ( ^ I (3 11) 

y diib Sijip + didjh + e ijk d k h J 

Keep in mind that the captical indices are lowered or raised by the Minkowski metric 
t]ab or its inverse. There are six scalar degrees of freedom altogether. The corresponding 
perturbed metric takes the form: 

/ 1 + 20 adi(w + w) \ 

y ^ \ adi{w + w) -a 2 ((1 - 2V)<% - 2didjh) )' K ' ' 

In general, w and w always affect the metric in terms of their combination w + w, which 
serves as a single degree of freedom, but as far as the vierbein is concerned, they are two 
independent degrees of freedom. Moreover, h does not present itself in the metric as well. 
To summarize, there are two degrees of freedom which do not contribute to the metric, 
and this is exactly what we ignore in the previous subsection. 

In longitudinal gauge, w = —w and h = 0. Since w has a mass dimension, we introduce 
a dimensionless quantity £ which is related to w by £ = aHw. As demonstrated in the 
appendix, one can obtain the perturbed equations up to first order as follows 

K 2 

E% : — 5p = a- 2 (l + f T )d 2 ^ - 12a~ 2 H 2 f TT d 2 C - 3H(1 + f T - l2H 2 f TT )^ + H</>), (3.13) 
E\: ^{p + P )d l v = -{l + f T ){d i ^ + Hd i (t>) + l2HHf TT d i ^, (3.14) 

K 2 

E\ - —(p + p)d i v = -{l + f T -12H 2 f TT ){dd + Hd i (j ) )-Aa- 2 Hf TT d i d 2 C, (3.15) 
Tr(^) : 

K 2 ( 1 

— 5p = (1 + fr) (iP + 3HiP + H<j> + 2H(t> + 3H 2 (j)) - -a~ 2 d 2 {^ - 0) 

+ f TT (-l2H 2 i> - 3QH(H + H 2 )i) - 12H 3 <]) - (60HH 2 + 36H 4 )(f) 

+a- 2 {8H + AH 2 )d 2 ( + 4a~ 2 Hd 2 C) 

+l2f TTT HH 2 (l2H(ip + Hcj)) - 4a" 2 9 2 c) , (3.16) 

E?j{i + j) : (1 + fr)di&(<l> - V) + 12HfTTdi&( = 0. (3.17) 

See Appendix ^ in detail. Note that in the above equations, the Parity-violating term 
h disappears, but w survives, even though w does not appear in the perturbed metric in 
longitudinal gauge. Compared to the first order equations in GR p8[|, we have an extra 
degree of freedom £ in the perturbed equations and one more equation is obtained. Our 
equations are self-consistent. On the other hand, for the trivial Lagrangian in which f(T) 
is a linear function of the torsion scalar T, all analysis should parallel to those in GR except 
for a re-scaled coupling constant, and the extra degree of freedom £ should disappear. Here 
we see that £ always appears in the company of Jtt- It is quite reasonable. 
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Since the matter Lagrangian is invariant under general coordinate transformation, the 
energy-momentum tensor should be conserved with respect to the Levi-Civita connection. 
Then one can find two equations which take exactly the same form as their counterparts 
in GR: 



bp + 3H(5p + Sp) + a~ 2 (p + p)d 2 v - 3(p + p)i> = 0, 
pdiv + {p + p)div + di5p +{p + p)di4> = 0. 



d_ 
~dt 



((p + p)div) + 3H(p + p)div + dip +(p + p)di4> = 0. 



(3.18) 
(3.19) 

(3.20) 



One can also derive the above two equations from Eq. (|3.13j - 3.17 ). From now on we will 
focus on the universe only filled with dust-like matter, namely p = Sp = 0, and all equations 



will be transformed to Fourier space. From Eq.( 3.19 ), one reaches a very useful relation 

v = -<f>. (3.21) 
Define a gauge invariant fractional matter perturbation 



where 



Pm 



6p m = 5p m - 3Hp m v 



(3.22) 



(3.23) 



is the gauge-invariant comoving matter density perturbation. It can also be interpreted 
as the density perturbation on spacelike hypersurfaces orthogonal to comoving worldlines. 
Considering Eq.(gjjp and Eq.(|3~T5|), we obtain 



Sp r , 



rV-(l + h). 



(3.24) 



From Eqs. ( |3.14|) , ( p. 17 ) and ( 3.22 ), the evolution of matter density perturbation is given 
by 

k 2 k 2 C 

6 m = ^v- YlRtifTT^^—. (3.25) 

a a ^Pra 

In order to get the evolution equation of 5 m , we need to work out the solution of £ as well. 

In this subsection we focus on the non-trivial case with f TT / 0. From Eqs. ( |3~H| ) 
and ( |3.15| ), we obtain 

k 2 



On the other hand, Eq. (3.17) can be written by 



= V - 12if-^-C- 
1 + /T 

Combing the above two equations, we find 

3(H - H 2 )ip - 3Hip 



c 



fc 2 

a? 



3mH 2 4¥r- 

l+JT 



(3.26) 



(3.27) 



(3.28) 
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In the subhorizon limit, £ ~ " if} <C V- One can expect that ( will play an important 
role on the evolution of perturbations at large scales. But here we only focus on the physics 
in the subhorizon limit and we have 

<p ~ ip. (3.29) 



It is the same as that in minimally coupled GR. Therefore, from Eq. ( 3.14 ), Hv ~ <p ~ 
ip 3> C and hence the second term on the right-hand side of Eq. fl3.25|) can be neglected. 
Taking the time derivative of Eq. ( |3.25| ) , one can obtain 

k 2 

5 m + 2H5 m + —0 = 0, (3.30) 
a z 



where Eq.( [3.2i ) is taken into account. Plugging Eq.( |3.29] ) into this equation and combining 



with Eq.( [3,24 ), the evolution equation of linear matter perturbation becomes 

d m + 2H5 m - 4:irG efi pm5m = 0, (3.31) 
where G e g is the effective Newton's constant which is related to G by 

GeS = T^T- (3-32) 

If f(T) is a linear function of T, namely f(T) = aT, the Newton's constant is just rescaled 
to be G e fi = G/(l + a), which is also constant in time. We can understand this result 
from the action ( p.7[ ) directly. In f(R) or scalar-tensor gravity, one can get an analogous 
equation for fractional matter perturbation 5 m with a redefinition of Newton's constant in 
the short wave-length limit p5fl . 

At the deep matter dominant era, if fx — 0, the solution of Eq. ( |3.3l| ) indicates that 
the matter density perturbation goes like 5 m oc a. Therefore 5p m = Pm^m ~ oT 2 and 
cj) ~ ip ~ constant in time. From Eq.( |3,26 ), C ~ a_1 which implies that ( decreases with 



respect to <f> and ifr. 

To summarize, we consider all the scalar degrees of freedom in this subsection and 
obtain the evolution equation of matter energy density perturbation in the subhorizon 
limit. In this limit, the extra degree of freedom £ is suppressed compared to 4> and tjj, but 
it removes the inconsistency in the former subsection. 

3.3 Numerical analysis 

For convenience, we define the growth as the ratio of the perturbation amplitude at some 
scale factor relative to some initial scale factor, D = 5 m {a) / $m{ a i) ■ The matter density 
perturbation 5 m is proportional to the scale factor a in the f(T) gravity with /y < 1 
during the matter era. We introduce a new variable g(a), namely 

9(a) (3.33) 
a 

which does not depend on a during the matter era; thus, the natural choice for the initial 



conditions are g(di) = 1, dg/d In a\ a=ai = 0. From Eq. (|3.3l|) , the equation for g{a) becomes 

d 2 q I , H \ do I H 4TvG e gp m , 

1 1 ' 1 • 1 " ' 1 g = 0. (3.34) 




dlna 2 \ H 2 dlna \ H 2 H 2 
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This equation reduces to that for dark energy scenario in GR with the same equation of 
state as the effective dark energy in f(T) if we replace G e flf by G. For a universe only filled 
with dust-like matter in f(T) gravity, we have 

H_ _ 3 1 + f/6H 2 + 2/ T 
H 2 

and 

4irG eS p m 3Q° m 



21 + / T -12tf 2 /tt' 



(3.35) 



-3 



H 2 - 2(1 + f T )H^ ' (3 ' 36) 

here the scale factor ao is normalized to be one and is the matter energy density 
parameter. The superscript '0' denotes that the variables are evaluated at present. 
For intance, we consider a power law model with 

f(T) = a{-T) n = a(6H 2 ) n , (3.37) 

where a can be determined by the present Hubble parameter and matter density parameter, 
namely 

a = (6H 2 ) 1 - n (l-n° m )/(2n-l). (3.38) 
The equation for the perturbation becomes 



d 2 9 
din a? 



+ 



+ 



4- 3 -± 
2 1 - 

3-21 
2 1 - 



h 



2n-2, 



1 



n° ) 



dg 

nh 2n - 2 (l-QP m )\ din a 

h 2 - 2 (i-nl) 



3nlh- 2 a- 3 



n 



h 2 " 



1 



2(1 



2n-l " / 



9 = 0, 



where h 



H/Hq which is governed by 

dh 2 -3h 2 + 3h 2n (l-n° m ) 



(3.39) 



(3.40) 



din a l-nh 2n - 2 (l-n^ r 
Since the above equations are too complicated to be solved analytically, we will use nu- 
merical method. In ||, Linder pointed out that this model can fit current observation only 
when n <C 1. Therefore we have 

'# n \2(l-») 



fr 



n 



i -<r. ) ( ^ 



(3.41) 



l-2n 

which is much smaller than one during the matter era because n <C 1 and Hq/H <C 1. 
Here we will adopt n = 0.1 and = 0.28 for numerical calculation. The initial moment 
should be taken during the matter era, e.g., Oj = 1/31 (i.e., z = 30). In addition, the initial 
condition of h(a) is h(a = 1) = 1. 

When / is a constant or n = 0, the term f(T) acts just as a cosmological constant. 
In this case, our result recovers that in ACDM in GR. For a universe filled with matter 
and dark energy whose equation of state is the same as that in f(T) gravity, the evolution 
equation of 5 m in the framework of GR can be obtained by replacing G e ff with G in ( 3.31] ) . 
Our numerical results are illustrated in Fig. |l[ Since fx > 0, the effective Newton's constant 
in f(T) gravity gets smaller than that in GR, and the gravitational interaction is weakened. 
That is why the over-dense perturbation in f(T) gravity grows slower than that in GR (See 
the blue dashed and black dotted lines). 
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z 

Figure 1: The numerical soluion of the evolution equation for the growth rate of matter pertur- 
bation in ACDM(the red solid line), power law f(T) model(the blue dashed line) and dark energy 
with the same equation of state in GR (the black dotted line). 

4. Discussions 

In this paper we derived the evolution equation for linear matter density perturbation in 
the framework of f(T) gravity and compared it to that in GR. We began our analysis 
from two aspects. One is based on a naive ansatz in which we chose the simplest vierbein 
and concluded that it leads to an inconsistency in Sec. 3.1. Though different vierbeins 
are related to each other by a Lorentz transformation, they may have different predictions 
because the Lorentz symmetry is broken down in f(T) theory. However, there is not a 
principle for us to choose a 'physical' vierbein. In Sec. 3.2, we proposed a strategy to solve 
this puzzle. We started with the most general perturbed vierbein and we found that an 
extra degree of freedom cure the inconsistency in the former case even though it does not 
appear in the perturbed metric. Finally, in Sec. 3.3, we figured out the growth factor in 
the power law /(T) model in detail and showed that the over-dense matter perturbation 
grows slower than that in GR due to the weakened gravity. 

When f(T) contains some nonlinear terms of T, f(T) gravity is not Lorentz invariant 
any more. However, the Lorentz symmetry should be preserved at least at small scales. 
This may require a stringent constraint on f(T) gravity. Once we take this constraint into 
account, whether f(T) theory can lead to an accelerated expansion of our universe is still 
an open question. We will come back to this question in the future. 
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A. Equations of motion for the perturbations 

In the appendix, we present some details of our calculations. From the definition in Eqs. 
(|3TT0|) and fu| ), one obtains the components of perturbed vierbein : 

e ° = 1 + 4>, ei° = adiw, 

erf = —d t w, e/ = a ^(1 — t/j)Sj — djd l h — ej m d n h^ , 
and its inverse e /J a'- 

e°o = 1 — 4>, e l o = a~ 1 d % w 1 

e°i = -8iW, e l j = a' 1 ((1 + ^)5) + dj&h + e 3 in d n h}j . 

Plugging the above equations into g^ u = r]ABG ll A e^ B , one gets the metric in Eq. ( 3.12| ). 
For simplicity, we choose the Newtonian gauge, where w = —w, and h = 0. In this case, 

T°oi = -di4> - adiw, 
T°ij = 0, 

T 0j = H8) - #j - e/ n d n h + a^d^w, 

r jk = d k ^ - d^5{ + e 3 m d k d n h - e k in d 3 d n h, 

and then 

S ° l = a" 2 aV, 

S ij = ~a- 2 e^ n d n ti, 

= -H5j + (ip + 2H(f> - \- l d 2 w)5{ + hr x di&w, 
Si jk = -^a- 2 e jkn did n h + ^a- 2 (d fc V - d k <f) - ad k w)5j - ^a~ 2 {d j ^ - - a&w)5 k . 
In addition, one can also calcualte perturbed torsion scalar, namely 

T = —QH 2 + 12H(tp + H(j>) — Aa~ l Hd 2 w. (A.l) 
Defining T a a = e p A^p a ^ one can easily obtain 

f\ = -p-5p + P <p, (A.2) 

f°i = a- 1 {p + p)d i v- pdi-w, (A.3) 

2*0 = -a~ 2 {p + p)d i v + aT^p&w, (A.4) 

= a - l (p + 8p)8) + a^pipSj + oT x pe 3 in dnh. (A.5) 
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Pluging these expressions into the equation of motion Eq. (|2.8| ), we find 
o: 

.2 



\p=\(T + f) + m 2 (l + f T ), (A.6) 



and 



!j(Sp + P<t>) = \{T + f)</> + (1 + /r)(3fty + 6# 2 - a" 2 ^) 

- 3H 2 f TT (l2H(iP + £T0) - 4a~ 1 #d 2 w). (A.7) 

Note that Eq. ( |A.6| ) describes evolution of homogeneous background and Eq.(A.7) equa- 



tion of motion for density perturbation 5 p. Combing Eq. (A.6) and Eq. ( |A.7| ) yields Eq. 



(3.13). 



K 2 1 

— (a- 1 (p + p)d l v-pd i w) = --{T + f)d i w + {l + f T ){-a- 1 d i ^-a- l Hdi(t ) -2,H 2 d l w) 

-1 0-2/0 „/. , rr« ,\ ^ -2rr2a.a2 



+ f T T{l2a- l H 2 {dii> + Hdi<p) - 4a- 2 H 2 did 2 w). (A.8) 



Considering Eq. ( |A.6| ), one can obtain Eq. ( |lg) . 



1 

y (-a~ 2 (p + p)9 t v + a- 1 p5 J u') = --a' 1 ^ + f)d l w - l2a~ 2 HHf TT (d l ip - aHd l w) 

+ (1 + / T ) (a" 2 ^ + #0) - a-^ij + 3F 2 )<9^Ya.9) 



-I(T + /)- (ii' + 3^ 2 )(l + / T ) + 12ij J H- 2 / TT , (A.10) 



2 



((^+p#*+pe/ n cU) 



1 



j(T + /)(#j + e/ n d n /i) + (1 + / T ) + 3HiP -Hifj- 3H 2 iP + H^> + 2H<j> + ZH 2 <f))5) 
-{H + 3H 2 )e/ n d n h - X -a- 2 d 2 {^ - <P)5) + ^a" 2 ^^ - 0)J 
+/tt ((-12H 2 ^ - 36H(H + # 2 )V> + l2HH 2 ip - l2H 3 <p - (60HH 2 + 36H 4 )cj))S : . 



3 

+a~ 1 (UHHd 2 w + 8H 3 d 2 w + 4i/ 2 <9 2 u;)^ + 12HH 2 ej in d n h - 6a~~ 1 HHdjd i w ' 
+12f T TTHH 2 (l2H(ip + H<t>) - ia~ 1 Hd 2 w) 5). (A.ll) 
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Plugging Eq. ( |A.10|) into ( |A.9|) , one reaches Eq. (|3.15|) . Using Eq. (|A.10|) , one can simplify 
Eq. flAlTD as 



(A.12) 



= (1 + f T ) Ui> + 3Hijj + H<j> + 2Hcf> + 3H 2 (f ) )5) - ^a- 2 d 2 (iP - 4>)5) + ^a^d^^ - 4>) 
+/rr ({-12H 2 ^ - 36H(H + H 2 )4> - l2H 3 <j) - (60HH 2 + 3QH 4 )(t))5) 
+ a - 1 (UHHd 2 w + 8H 3 d 2 w + m 2 d 2 w)5) - Qa^HHd^w^j 
+12 f TTT HH 2 (l2H(ip + H(j)) - Aa- 1 Hd 2 w^j 8). 

Taking the trace of Eq.( |A.12 ), one obtains 

k 2 { 1 

— Sp = (1 + / T ) ( (V + 3HiP + + 2H<P + 3H 2 et>) - -a- 2 d 2 {^ - <j>) 

+/tt (-12H 2 i) - 3QH(H + H 2 )ip - l2H 3 <p - (60HH 2 + 36H 4 
+a~ 1 {12HHd 2 w + 8H 3 d 2 w + AH 2 d 2 w)J 
+12f TTT HH 2 (l2H(ip + H(j>) - Aa' 1 Hd 2 w^j . 



(A.13) 



Noticing that <^ = aHw, we can rewrite this equation in terms of Eq.( 3.16 ). In addition, by 
combining Eq.( [A.12 ) and Eq.(A.13), one can obtain Eq.( 3.1?1 ) after expressing w in terms 
of the dimensionless quantity Here we want to stress that the terms with h in ( |A.ll ) 
are cancelled and thus h does not show up in the perturbation equations in Sec. 3.2. 
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